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Abstract 

■ A steep maximum occurs in the Wroblewski ratio between strange and 
<^ ■ non-strange quarks created in central nucleus-nucleus collisions, of about 

A=200, at the lower SPS energy yfs ~ 7 GeV. By analyzing hadronic mul- 
tiplicities within the grand canonical statistical hadronization model this 

■ maximum is shown to occur at a baryochemical potential of about 450 
MeV. In comparison, recent QCD lattice calculations at finite baryochem- 
ical potential suggest a steep maximum of the light quark susceptibility, to 
occur at similar hb, indicative of "critical fluctuation" expected to occur 
at or near the QCD critical endpoint. This endpoint hat not been firmly 

I/-) | pinned down but should occur in the 300 MeV < fi B < 700 MeV inter- 

CN ■ val. It is argued that central collisions within the low SPS energy range 

should exhibit a turning point between compression/heating, and expan- 

Q sion/cooling at energy density, temperature and (jlb close to the suspected 

critical point. Whereas from top SPS to RHIC energy the primordial dy- 
namics create a turning point far above in e and T, and far below in fig. 

^ ' And at lower AGS energies the dynamical trajectory stays below the phase 

boundary. Thus, the observed sharp strangeness maximum might coincide 
with the critical at which the dynamics settles at, or near the QCD 
r*| , endpoint. 
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Bulk hadron production systematics in central nucleus-nucleus collisions at 
relativistic energy is, overall, well reproduced by a statistical Hagedorn hadronic 
freeze-out model. A grand canonical version of this model captures the various 
hadronic species multiplicities, per collision event, from pions to omega hyper- 
ons, in terms of a few universal parameters that describe the dynamical stage in 
which the emerging hadronic matter decays to a quasi-classical gas of free reso- 
nances and hadrons [1, 2, 3]. The grand canonical parameters are temperature 
T, volume V and chemical potential \x. They capture a snapshot of the fireball 
expansion within the narrow time interval surrounding hadronic chemical freeze- 
out, which thus appears to populate the hadron/resonance mass and quantum 
number spectrum, predominantly, by phase space weight [4, 5] thus creating an 
apparent thermal equilibrium state prevailing in the produced hadron-resonance- 
population. This chemical equilibrium instantaneously decouples from fireball 
expansion surviving further (near isentropic) processes. It can thus be retrieved 
from the finally observed hadronic multiplicities, by state of the art grand canon- 
ical model analysis. This analysis succeeds from AGS, via SPS, to RHIC energy 
[1, 2, 3, 5]. 
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Statistical model analysis is also applicable to elementary collisions, p + p, 
p+p, and e + e~ annihilation as was shown by Hagedorn [6] and, more recently, by 
Becattini and Collaborators [7, 8]. The canonical version of ensemble analysis is 
applicable here. Mutatis mutandis the same hadrochemical equilibrium feature is 
being attested, emphasizing the statement that the apparent equilibrium does not 
arise from a thermodynamical inelastic rescattering cascade toward equilibrium - 
there is essentially none in elementary processes - but should stem directly from 
the QCD hadronization process occuring under phase space dominance [4, 5]. 

The crucial difference between elementary and central nucleus-nucleus resides, 
in statistical model view, in a transition from canonical to grand canonical or- 
der in the ensuing decoupled hadronic state. This transition was studied by 
Cleymans, Tounsi, Redlich et al. [9]. Its main feature is strangeness enhance- 
ment. Comparing the strange to non-strange hadron multiplicities in elementary, 
and in central nucleus-nucleus collisions at similar energy, one observes an in- 
crease of the singly strange hyperons and mesons, relative to pions, of about 2-4, 
and corresponding higher relative enhancements of multiply strange hyperons 
[10, 11, 12], ranging up to order-of-magnitude enhancement. In the terminology 
of Hagedorn statistical models, strangeness is suppressed in the small system, 
canonical case, of elementary collisions (due to the dictate of local strangeness 
conservation in a small "fireball" volume), whereas it approaches flavour equipar- 
tition in large fireballs due to the occurence of quantum number conservation, on 
average only, over a large volume - as reflected by the global chemical potential 
featured by the grand canonical ensemble: "strangeness enhancement" occurs as 
the fading-away of canonical constraints. 

From statistical model analysis we obtain a more general view of strangeness 
relative to non-strangeness production than is provided by considering individual 
strange to non-strange production ratios, like K/ir, Q/ir etc., from p+p to central 
A+A. The model quantifies strange to non-strange hadron/resonance production 
by means of Wroblewski quark counting at hadronic freeze-out [13]. It determines 
the so-called Wroblewski-ratio, 

\ s = 2(<s> + <s>) 

< u > + < u > + < d > + < d > 

which quantifies the overall strangeness to non-strangeness ratio at hadronic 
freeze-out. Strangeness enhancement (i.e. removal of strangeness suppression in 
elementary collisions) is quantified, by such an analysis, to proceed from \ s ^0.25 
in elementary collisions, to A ~ 0.45 in central nucleus-nucleus collisions [2, 3]. 

Now we turn to the point of the present study. From a recent energy scan 
conducted at the SPS by NA49, studying hadron multiplicities from ^=7 to 17 
GeV, a steep maximum was observed [14] in the K + /n and A/n ratios in central 
Pb+Pb collisions, as shown in Figs.l and 2. As the K + and A channels carry 
most of the total < s > + < s > content, this experimental result indicates a 
kind of "singularity" in the strange to non-strange production ratio, from AGS 
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Figure 1: Energy dependence of < K + > / < 7r + > ratio for central Pb+Pb 
(Au+Au) collisions (upper points) and p+p interactions (lower points). 

to RHIC energy. It appears to be unlikely that state of the art hadronic or 
partonic quasi-classical microscopic transport models should exhibit such non- 
smooth behaviour which is also absent in p + p collisions. In order to generalize 
the new NA49 data, away from consideration of individual channel strange to 
non-strange multiplicities, Becattini et al. [3] analyzed the yfs dependence of the 
Wroblewski-parameter X s in the grand canonical statistical hadronization model. 
Their result is shown in Fig. 3 which gives X s as a function of the chemical potential 
fi B - From top AGS energy (at p, B ~ 550 MeV) to RHIC energy (p, B < 50 MeV) 
one perceives an average A s of about 0.45 ± 0.08 whereas a steep excursion is seen, 
to A s = 0.6 ± 0.1, at /ie=440 MeV. This point corresponds to the steep maxima 
observed in Figs.l and 2, to occur at SPS fixed target energy of 30 GeV/A in 
central Pb+Pb collisions, corresponding to y/s=7. 3 GeV. The steep singularity 
of the NA49 K + /it and A/tt ratios at this y/s thus reflects in a maximum of the 
Wroblewski A s , derived from grand canonical analysis. For the sake of clarity 
we state here that the latter analysis just takes note, in a snapshot at hadronic 
freeze-out (to a decoherent hadron/resonance gas) of the general multiplicity 
order prevailing at the instant of hadronization. The observed X s maximum 
should, therefore, present a hint that hadronization at y/s ~ 7 GeV should occur 
under influences, absent at energies above and below. Moreover, NA49 has shown 
recently [15] that the event-by-event fluctuation of the ratio (K + + K~)/(tt + + 
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Figure 2: Energy dependence of the <A>/<7r>(<7r >= 1.5 (< n~ > 
+ < 7r + >)) ratio for central Pb+Pb (Au+Au) collisions (upper points) and p+p 
interactions (lower points). 



ii ) measured in central Pb+Pb collisions increases steeply toward y/s = 7 GeV 
whereas it was formerly found [16] to amount to be below 4%, at top SPS energy, 
y/s = 17.3 GeV. 

We thus propose that the dynamical trajectory of central Pb+Pb collisions 
comes close to the critical point of QCD, at or near y/s = 7 GeV. This point 
has been expected to occur on the line in the T, \x& plane which describes the 
boundary between the hadronic and partonic QCD phases [17]. Along that line 
the phase transition is expected to be a crossover at \ib < I^b, become second 
order at fis = \f B i and first order for \ib > H B - At [i B = ^ B and T = T c we 
thus expect phenomena analogous to critical opalescence. Recent QCD lattice 
calculations succeeded in an extrapolation to finite chemical potential [18, 19], 
thus making a first prediction for the phase boundary line and, in particular, the 
critical point - albeit with considerable uncertainty as was discussed by Redlich 
at QM04 [20]. This uncertainty stems, firstly, from the uncertainty in the ex- 
trapolation to finite Hb but, secondly, from the unphysical (high) strange quark 
mass employed in these lattice calculations which, at present, place the critical 
point somewhere in the interval 500 MeV < fi c B < 700 MeV. Redlich argued that 
it should move to considerably lower ^ once the s-mass can be chosen closer 
to the physical quark mass [20]. This expectation was substantiated by recent 
lattice calculations which show that the critical point might move downward in 
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Figure 3: Dependence of the X s parameter on baryochemical potential extracted 
from the fits to hadron multiplicities in central Pb+Pb (Au+Au) collisions at 
AGS, SPS and RHIC energies. The lines show the dependence expected for 
different values of the 7 S parameter. 



fiB once more realistic quark masses are employed [21]. From Fig. 3 we see that 
the strangeness maximum at = 7 GeV corresponds to fis ~ 440 MeV and 
thus quite close to the expected \i c B position. Furthermore the energy density at 
the phase boundary is estimated by lattice QCD to be rather low [22] (e < 1 
GeV/fm 3 ). 

Central collisions of heavy nuclei at SPS energy exhibit a general cycle of 
initial compression and heating which is followed by a maximum energy density 
stage which then turns into expansion and cooling [23]. The quantities charac- 
terizing the overall system dynamics, such as volume and energy-entropy density 
etc. change very rapidly except during the high density stage which acts analo- 
gous to a classical turning point. Also it is only during this stage that relaxation 
times can be of comparable magnitude to the system evolution time scales (like 
volume-doubling etc.). One can thus begin only here treating the system dynam- 
ical evolution in terms of energy density, temperature and chemical potential - 
thus defining a dynamical trajectory in the T, ^b plane. Only if this turning 
point coincides closely with the QCD critical endpoint one could expect to ob- 
serve substantial critical phenomena. Now it is well known that the maximum 
energy density in central collisions of mass 200 nuclei amounts (Bjorken estimate) 
to above 2 GeV/fm 3 at top SPS [24], and to about 5 GeV/fm 3 at RHIC [25] ener- 
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gies, thus overshooting, by far, the critical QCD energy density. The system will 
thus cross the phase coexistence line, upon re-expansion, whilst already undergo- 
ing rapid expansion. Furthermore, the chemical potential is certainly well below 
300 MeV at the time of hadronization.The evolution will thus miss the critical 
point at top SPS, and RHIC energies; and at much lower AGS energies the dy- 
namics falls into the \xb >500 MeV domain but the energy might not suffice to 
reach the phase boundary. In summary we may indeed expect that the dynami- 
cal evolution reaches its energy density plateau phase near the expected critical 
point (i.e. at energy density just below or at 1 GeV/fm 3 , and at hb between 300 
and 500 MeV) somewhere in between maximum AGS and minimum SPS energy. 
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Figure 4: The quark number susceptibility calculated within lattice QCD as 
function of temperature (relative to transition temperature) for different values 
of quark chemical potential. 



We now turn to the final point of our line of argument by returning to the 
lattice results [19, 20] at finite They see a steep maximum of the quark 
number susceptibility 

Xu ' d = T2{ d(fi/ry (2) 

occuring at T=T C = 150 MeV and /is = 3/iQ=3 T = 450 MeV. We reproduce 
the Bielefeld- Swansea results in Fig. 4 which also shows the calculations for fig = 
225 MeV and hb = (essentially corresponding to top SPS and RHIC energies, 
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respectively). The latter exhibit no susceptibility peak but a smooth transition 
from T < T c to T > T c . As Xu,d can also be written as 

_ T 2 ( 6 , n u + n d 

Xq ~ 1 <W) + KmtY [6) 

we see that the peak in the susceptibility implies a maximum fluctuation of the 
quark number densities n u and n^. We interpret this result as an indication 
of critical fluctuation occuring in the vicinity of the critical endpoint implicitly 
present in this calculation. Directly at fi c the susceptibility would diverge. The 
critical point in this calculation must thus be near /ib — 450 MeV and T = 150 
MeV. This, in turn, is very close to the parameters of the grand canonical model 
at the strangeness maximum, y/s = 7 GeV (Fig. 3). 

As to the relation between the susceptibility maximum of lattice QCD and the 
strangeness maximum observed by NA49 (at which the Wroblewski parameter X s 
exhibits an anomaly), Gavai and Gupta [26] have suggested the relationship 

A, = (4) 

X.u ~r Xd 

which appears to offer a direct link. In fact they obtain X s = 0.48 from a lattice 
calculation at zero /if,: closely coinciding with the value observed at top SPS and 
RHIC energy (Fig. 3). Unfortunately, though, their result refers to Hb = 0, and 
the Bielefeld-Swansea calculations at finite [19] arc in two-flavour QCD only. 
A prediction for Xs at hb ~450 MeV, or, more generally, a full three-flavour 
lattice treatment of the vicinity of the critical point is required to finally assess 
the line of argument of the present note. If proven correct we would encounter, 
here, the first direct reflection of QCD in nucleus-nucleus collision data. 

The author acknowledges stimulative discussions with F. Becattini, R. V. 
Gavai, F. Karsch and K. Redlich. 



References 

[1] F. Becattini, M. Gazdzicki and J. Sollfrank, Nucl. Phys. A638 (1998) 403; 
J. Cleymans and K. Redlich, Phys. Rev. Lett. 81 (1998) 5284. 

[2] P. Braun-Munzinger, I. Heppe and J. Stachel, Phys. Lett. B465 (1999) 15; 
P. Braun-Munzinger, D. Magestro, K. Redlich and J. Stachel, Phys. Lett 
B518 (2001) 41; 

P. Braun-Munzinger, K. Redlich and J. Stachel, nucl-th/0304013; 

P. Braun-Munzinger, J. Cleymans, H. Oeschler and K. Redlich, Nucl. Phys. 

A697 (2002) 902. 

[3] F. Becattini, M. Gazdzicki, A. Keraenen, J. Manninen and R. Stock, Phys. 
Rev. C69 (2004) 024905. 



7 



R. Hagedorn, Nucl. Phys. B24 (1979) 93; 

J. Ellis and K. Geiger, Phys. Rev. D54 (1996) 1967. 

R. Stock, Phys. Lett B456 (1999) 277; 
R. Stock, hep-ph/0312039. 

R. Hagedorn, Nuovo Ciemento 35 (1965) 395. 

F. Becattini, Z. Phys. C69 (1996) 485. 

F. Becattini and L. Ferroni, hep-ph/0307061. 

J. Cleymans et al., Phys. Rev. C56 (1997) 2747; 

A. Tounsi and K. Redlich, J. of Physics G28 (2002) 2095. 

S. V. Afanasiev et al., NA49 Coll., Phys. Rev. C66 (2002) 054902; 
A. Mischke et al., NA49 Coll., Nuc. Phys. A715 (2003) 453. 

L. Ahle et al., E802 Coll., Phys. Rev. C60 (1999) 064901. 

F. Antinori et al., NA57 Coll., Nucl. Phys. A698 (2002) 118; 
F. Antinori, talk at QM04. 

A. K. Wroblewski, Acta Phys. Pol. B16 (1985) 379. 

V. Friese et al., NA49 Coll., nucl-ex/0305017; 
M. Gazdzicki, talk at QM04. 

Ch. Roland et al, NA49 Coll., talk at QM04. 

S. V. Afanasiev et al., NA49 Coll., Phys. Rev. Lett. 86 (2001) 1965. 

M. G. Alford, K. Rajagopal and F. Wilczek, Phy. Lett B442 (1998) 247; 
Y. H. Atta and T. Ikeda, Phys. Rev. D67 (2003) 014028. 

Z. Fodor and S. D. Katz JHEP 0203 (2002) 14. 

R. V. Gavai and S. Gupta, hep-lat/0303013; 
C. R. Allton et al., hep-lat/0305007. 

K. Redlich, talk at QM04. 

Z. Fodor and S. D. Katz, hep-lat/0402006; 

F. Karsch et al., Bielefeld-Swansea lattice Collaboration, hep-lat/0309116. 
F. Karsch, Nucl. Phys. A698 (2002) 199. 
R. Stock, Phys. Reports, 135 (1986] 259. 



8 



[24] T. Alber et al, NA49 Coll., Phys. Rev. Lett 75 (1995) 3814. 
[25] G. Roland, talk at QM04. 

[26] R. V. Gavai and S. Gupta, Phys. Rev. D (2002); 
R. V. Gavai, talk at QM04; 
R. V. Gavai and S. Gupta, hep-lat/0303013. 



9 



